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We study controlled phasegates for ultracold atoms in an optical potential. A shaped laser pulse 
drives transitions between the ground and electronically excited states where the atoms are subject 
to a long-range 1/R 3 interaction. We fully account for this interaction and use optimal control 
theory to calculate the pulse shapes. This allows us to determine the minimum pulse duration, 
respectively, gate time T that is required to obtain high fidelity. We accurately analyze the speed 
limiting factors, and we find the gate time to be limited either by the interaction strength in the 
excited state or by the ground state vibrational motion in the trap. The latter needs to be resolved 
by the pulses in order to fully restore the motional state of the atoms at the end of the gate. 
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I. INTRODUCTION 



The physical realization of a quantum computer re- 
quires the implementation of a set of universal gates [l[ . 
The most difficult part is generally the two-qubit gate 
since it involves interaction between two otherwise iso- 
lated quantum systems. In proposals for quantum com- 
puting with ultracold neutral-atom collisions @, Q , the 
two-qubit gate operation involves atomic motional de- 
grees of freedom [j,[5j], most often following adiabatic pro- 
cesses. This implies frequencies much lower than those 
characteristic of the trap, typically around a few tens 
of kHz. When long-range interactions, like dipole-dipole 
forces between Rydberg atoms [HI], are employed, the 
relevant energy scales are larger and gate speed can in 
principle reach a few GHz. 

Here, we study the limits to the two-qubit gate opera- 
tion time for resonant excitation of two ultracold atoms 
into an electronically excited molecular state. This may 
be a low-lying state, like those used for photoassociating 
two atoms to form a molecule d G3, or high-lying Ryd- 
berg states In this scenario, the system dynamics 
becomes more complex, involving motion under the in- 
fluence of the excited state potential. A high-fidelity gate 
can then no longer be designed "by hand" . Fortunately, 
since any gate operation corresponds to a unitary trans- 
formation on the qubit basis, its implementation can be 
formulated as a coherent control problem [IJ [l2| ■ Solu- 
tions to the control problem can be found theoretically 
within the framework of optimal control [l3r{T5| . Such an 
approach has been explored theoretically for molecular 
quantum c omp uting with qubits encoded in vibrational 
states [111 . Ilq . Experimentally the control problem is 
solved using femtosecond laser pulse shaping combined 
with feedback loops [ijj . Implementation of single-qubit 



gates with shaped picosecond pulses has recently been 
demonstrated for a qubit encoded in hyperfine levels of 
an atomic ion (l8| . In this experiment, the fundamental 
limit to the gate operation time is set by the inverse of 
the hyperfine splitting since the hyperfine dynamics are 
required to realize arbitrary qubit rotations. The gate 
duration is in particular much shorter than the period of 
atomic motion in the trap (l8j . 

Our goal is to implement a controlled phasegate, 



6 = diag (e* 1,1,1) 



(1) 
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between two qubits carried by neutral atoms. Shaped 
short laser pulses are used to drive transitions into an 
electronic state where the two atoms are interacting. Ob- 
viously, the minimum gate operation time will depend 
on the interaction strength. A second timescale comes 
into play because the interaction couples electronic and 
nuclear dynamics, inducing vibrational excitations of the 
two atoms. The vibrational period of the trap might thus 
also affect the minimum gate operation time. 

The limit to how fast a quantum gate can be per- 
formed is closely related to the minimum time it takes 
for a quantum system to evolve from an initial state to 
an orthogonal state, see for example [19j and references 
therein. This bound has been named the quantum speed 
limit; it is given in terms of the average energy and en- 
ergy uncertainty of the quantum evolution [201 ] - Since 
evaluation of the bound requires knowledge of the full 
spectrum, one can typically evaluate it analytically only 
for simple model systems. Numerically, the bound can be 
determined using optimal control where the breakdown 
of convergence indicates that the quantum speed limit 
has been reached [2l|. The following procedure allows 
us to determine the quantum speed limit for our desired 
two-qubit gate: We first set the gate operation time to a 
sufficiently large value to obtain a high fidelity implemen- 
tation by optimal control. Then we reduce the gate time 
until the optimization algorithm cannot anymore find a 
high-fidelity solution. Analysis of the fairly complex sys- 
tem dynamics reveals which part of the overall system 
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FIG. 1: (color online) Two calcium atoms in neighboring sites 
of an optical lattice (a), separated by distance d. The lattice 
potential can be taken in a harmonic approximation, allowing 
to separate the motion of both atoms into the center-of-mass 
motion and the interatomic motion (b). A phasegate is im- 
plemented by applying a shaped short laser pulse. 



dynamics limits the gate operation time. 

We consider neutral atoms trapped in optical tweezers 
or in neighboring sites of an optical lattice (Fig. [IJ . The 
sites to the right and left of the atom pair are assumed to 
be empty. We first consider alkaline-earth atoms which 
possess extremely long-lived excited states, cf. Ref. [22[ 
and references therein. The qubits can therefore be en- 
coded directly in the electronic states, that is, the atomic 
ground state 1 S$ together with the 3 Pi clock-transition 
state form a qubit. The laser is slightly detuned from the 
dipole allowed atomic transition between the 1 Sq ground 
state and the 1 Pq excited state, exciting the atom pair 
into the i? 1 molecular state. When both atoms are 
in the ground state, their interaction is of van-der-Waals 
type and practically zero at the trap distance. Due to dif- 
ferent exchange interaction in the electronically excited 
state, the B 1 ^ state scales as 1/R 3 at long range. This 
interaction may be employed to entangle the qubits pro- 
vided the time that the atom pair resides in this state is 
much shorter than its lifetime of a few nanoseconds. The 
interaction strength is determined by the distance of the 
atoms. For realistic lattice parameters, d > 200 nm, the 
interaction is too weak to generate entanglement in a suf- 
ficiently short time. On the other hand, we would like to 
probe what are the factors that limit the achievable speed 
on a general level. To this aim, we first explore a regime 
that cannot be realized in experiments, but presents a 
clear separation of timescales that allows for a clearer 
interpretation of the dynamics. We start assuming a fic- 
titious distance of d — 5 nm and a correspondingly un- 
realistic atomic trap frequency of 400 MHz. This allows 
us to identify the limiting factors for fidelity and gate 
time. Formally, our Hamiltonian is equivalent to the one 
yielding a Rydberg phasegate for alkali atoms @. We 
therefore study in a second step the implementation of 
a phasegate based on very strong dipole-dipole interac- 
tion in the excited state for realistic lattice spacings. In 
particular, we seek to answer the question whether ex- 
citation of the vibrational motion can be avoided if the 
excited state interaction is strong enough to allow for 
very short gate pulses. 



FIG. 2: Levels and laser-induced transitions for two qubits 
encoded in electronic states of two calcium atoms. 
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FIG. 3: (color online) Potential energy curves describing the 
interaction of two calcium atoms in the five lowest electronic 
states. The asymptotic values corresponding to R — oo are 
indicated on the right. 



The paper is organized as follows. Section [XT] intro- 
duces our model for the two atoms and summarizes how 
quantum gates can be implemented using optimal control 
theory. The numerical results are presented in Sec. IIIII 
with Sec. IIII Al devoted to generation of entanglement for 
two calcium atoms via interaction in the B 1 ^^ state. A 
generic dipole-dipole interaction, —C3/R 3 , is considered 
in Sec. IIII 51 where we determine the gate operation time 
for varying C3. We draw our conclusions in Sec. IIVI 



II. THEORETICAL APPROACH 

A. Modeling Two Atoms in an Optical Lattice 

We consider the following qubit encoding in a single 
calcium atom: The So ground state corresponds to the 
qubit state |0) and is used to define the zero of energy. 
The 3 Pi first excited state, taken to be the qubit state 
|1), then occurs at E\ = 15210cm" 1 . We consider the 1 P\ 
level as auxiliary state, |o) , with energy E a = 23652 cm" 1 . 
For two atoms we then obtain nine electronic states |00), 
|01), \0a), |10), |11), I la), \a0), |al), with asymptotic en- 
ergies Eij = Ei + Ej for the states \ij) — 0,1, a). 
This is depicted in Fig. [5] With the laser tuned close to 
the transition |0) <H- \a), no further electronic states will 
be resonantly populated. The potentials describing the 
interaction between the two atoms in the five lowest elec- 
tronic states are shown in Fig. [3) The AT 1 E+ ground state 
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potential corresponding to the two-qubit state |00) shows 
a 1/R 6 behavior at long range such that the atoms are 
effectively non-interacting at any relevant distance, while 
the auxiliary B 1 ^^ state corresponding to \0a) and |a0) 
goes as 1/R 3 . The remaining potentials corresponding to 
the two-qubit states |oo), |al), |la), |11), |10), and |01), 
are also essentially zero at the relevant distances. The 
potentials and transition dipole moment functions em- 
ployed in the following calculations are gathered from 
Refs. [H, l24j . For a Rydberg phasegate with alkali 
atoms, the qubit encoding is the standard one in hyper- 
fine levels of the electronic ground state; and the auxiliary 
level corresponds to the Rydberg state (we omit in this 
description any intermediate level that might be needed 
for near-resonant two-photon excitation of the Rydberg 
state). In this case, the ground state potential is sim- 
ply set to zero, while the auxiliary state is modeled by a 
generic C3/R 3 potential. 

The parameters of the optical lattice are chosen such 
that the motion of the atoms is restricted to one spa- 
tial dimension 3> u)q). Approximating the resulting 
trapping potential for the two atoms by two displaced 
harmonic oscillators, we can separate the center of mass 
coordinate which is subsequently integrated out, and the 
relative coordinate, i.e. the interatomic distance R. 

The Hamiltonian of our model, taking into account 
the nine two-qubit states, the motion in the trap, and 
the interaction between the atoms, is given by 

H = H lq ® l lq ® Ifl + llg ® Hlq ® Ifi + Hint 

= \ ik )( ik \ [ f + Vtrap(i?) + VBoiR) + Eik 

i,k 

+e(t) £ WW + \ki)(kj\] ® A«(i2). (2) 

Here, fii 9 denotes the Hamiltonian for a single three- level 
atom, li q is the identity in SU(Z), and 1r represents the 
identity for the motional degree of freedom. The inter- 
action Hamiltonian H; nt contains the Born-Oppenheimer 
potentials \/bo, and the laser pulse e(t) couples to the 
transition dipoles p,^. The indexes i,j, k each run over 

0. 1, a. The Hamiltonian, Eq. ([3]), is represented on a 
Fourier grid with variable step size [25l427j | . 

A two-qubit gate is successfully implemented if the four 
basis states |00), |01), |10), |11) transform according to 
the desired unitary transformation 0. Initially, the wave- 
function for the motional degree of freedom is given in 
terms of two atoms in the ground state of the displaced 
harmonic oscillators, ifo(R) — (R\ipo). Hence we con- 
sider four initial states given by \ij(R)) = \ij) Cg> |vo)> 

1, j = 0,1. The dynamics induced by a laser pulse will 
populate other states in the full 3 x 3 x A^-dimcnsional 
Hilbcrt space (with Nr the number of grid points to rep- 
resent the motional degree of freedom); that is, it will 
induce internuclear motion leading out of the logical sub- 
space. In order to calculate the gate operation in the log- 
ical subspace at the end of the pulse, the final states are 



reduced to the logical basis by tracing out the motional 
degree of freedom after projection onto \tpo). 



B. Optimal Control Theory for a Two-Qubit Gate 

In order to implement the target operation defined 
in Eq. ([1]), a suitable pulse e(i) must be found that 
drives the system evolution such that U(T, 0;e) = 
0. Optimal control treats the fidelity F, F e [0,1], 
which measures how close the evolved two-qubit basis 
states U(T, 0; e) \ij(R)) come to the desired target states 
O \ij(R)), as a functional of the control field e. Allow- 
ing for additional constraints such as minimization of the 
integrated pulse energy leads to the total functional J, 



J = —F 



g(e)dt, 



(3) 



which is to be minimized. Variation of the functional J 
with respect to the evolving two-qubit basis states and 
the control field yields a set of coupled optimization equa- 
tions that need to be solved iteratively jl2Hl5ll. We use 
the linear Krotov algorithm as outlined in Ref. [12[ to per- 
form the optimization. Starting with an arbitrary guess 
pulse e( k \t), k = 0, the algorithm sequentially updates 
the pulse at every point in time to yield an optimized 
pulse e( k+1 \t) that is guaranteed to improve the target 
functional J in each step of the iteration k. We take the 
running cost g(e) in Eq. ^ to be the change in integrated 
pulse energy, 



9(e) 



S(t) 



(*) 



Jk) 



(*) 



(4) 



rather than the integrated pulse energy itself. This guar- 
antees that, as we approach the optimum, 17(e) goes to 
zero such that the minimum of J becomes equal to the 
maximum of F [12]. In Eq. (j4]), a denotes an arbitrary 
positive scaling parameter, and the shape function S(t), 



S(t) = sin 2 (7rf/T) 



(5) 



ensures that the pulse is switched on and off smoothly. 

A possible starting point to define the fidelity is given 
by the complex- valued inner matrix product [12J, 



6 f u 



ij(R) 



(6) 



Out of the several choices for obtaining a real-valued 
functional, we employ 



F = [r] , 



(7) 



which is sensitive to a global phase [12]. The Hilbert 
space of the optimization problem can then be reduced 
to the 8iV#-dimensional subspace since the 1 11) level has 
no coupling to any other level, see Fig. [5] The evolution 
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of the ll(i?)) state cannot be controlled by the laser 
pulse but it is known to be 

U(T,0;e) \ll(R)) = e^ T \ll(R)) ; T = E n T/H. 

Including the information about the phase 0t explicitly 
as a global phase for all target states O \ij(R)), the evolu- 
tion of |ll(i?)) can be omitted from the system dynamics. 
Since entanglement is generated only by the dynamics 
of the states in the left-most column of Fig. [21 and the 
remaining levels remain uncoupled from the vibrational 
dynamics, the system dynamics can be further separated 
into those of the four molecular states in the left-most col- 



umn of Fig. [2] and those of a two-level system represent- 
ing the states in the middle columns of Fig. [21 Care must 
then be taken to extract the true two-qubit phase from 
the evolution of the |00(i?)} state which contains two- 
qubit and single-qubit contributions. This is described in 
Appendix [A] The numerical results presented below are 
obtained both within the 8./Vfj-dimensional model and 
the 4iVR-dimensional model plus two-level system, and 
the optimal pulses have been cross-checked. 

For the fidelity defined in Eq. ((TJ) and the running cost 
given in Eq. ([4]), the optimization equations read [12| 



Ae(f) 



e (*+D( t )_ e (*)(t) = ^*)am 
2a 

(ij(R)\6 f &(t,T;e^) ) 



= U(t,0;e 



MR)) 



(^(R^)\^(R)\H W ( R ^ 

«,i=o,i 



(8) 

(9) 
(10) 
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Here, (^f^j(R;t)\ denotes the backward propagated tar- 
get state \ij(R)) at time t. The backward propagation 
is carried out using the old field, e^(t). Uf{™{R;t) 

represents the forward propagated initial state \ij(R)) at 
time t. The new field, e^ k+1 \t), is employed in the for- 
ward propagation. (R;t)\ and y{™(R;t)\ are ob- 
tained by solving the time-dependent Schrodinger equa- 
tion numerically with the Chebychev propagator [28]. 
The time is discretized in n t steps of width At, between 
and T. Since no rotating wave approximation is em- 
ployed, At has to be fairly small (0.025 fs to 0.05 fs). 

For the desired gate implementation, there are two as- 
pects to the optimization problem: On one hand, the 
two-qubit phase x has to be realized. This is possible due 
to the interaction between the two qubits in the electron- 
ically excited auxiliary state. On the other hand, control 
over the motional degree of freedom has to be exerted. 
That is, at the end of the gate operation, the motional 
states \ij(R)) = \ij) ® \<po) have to be fully restored, 
except for the phases 0jj. Final wavefunctions contain- 
ing contributions from eigenstates other than \ij) (g) \ipo) 
imply leakage from the quantum register. These two as- 
pects of the optimization result can be quantified inde- 
pendently, and allow for a more thorough analysis of the 
solutions to the control problem than just the fidelity. 

The success of control over the motional degree of free- 
dom for the ground state is measured by projecting the 
final state, U(T, 0; e opt ) |00(i?)), onto the desired state, 



oo 



00{R) G(T,0;e opt ) 00(A) 



(11) 



The phase acquired by each of the propagated basis states 
is given by 

<fa = arg ((ij(R) |0(T,0;e**)| ij(R))) . (12) 

The phases 4>ij contain both single-qubit and two-qubit 
contributions. This is due to the weak molecular in- 
teraction which corresponds to small detunings of the 
laser from the atomic transition line. The pulse there- 
fore drives single-qubit purely atomic local transitions in 
addition to true two-qubit nonlocal transitions into the 
molecular state. The Cartan decomposition of a two- 
qubit unitary into local and non-local contributions pro- 
vides a tool to extract the desired non-local phase x from 
the <pij . It states that any two-qubit gate can be decom- 
posed as U = ki Ak2, where ki and l<2 are local operations 
acting on the single qubits, and A is a non-local opera- 
tor. The non-local contribution is characterized by local 
invariants [29]. For the phasegate considered here, the 
local invariants are calculated according to Ref. [29| to 
be 



Gl (0) 
G 2 (0) 



= cos 



»00 



001 - 010 + 4>ll) 



2 + cos [0i 



00 



'01 



'10 



On the other hand, for the desired target gate O defined 
in Eq. (TTJ), the local invariants are given by 



Gi 



( 6 ) 



cos 



G 2 6 = 2 + cos[ X ] 



We can thus identify the non-local phase, 
X = 0oo - 0oi - 0io + 0ii ■ 



(13) 
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The non-local character of the implemented gate U can 
also be measured in terms of the entangling power or con- 
currence [13] ■ For the controlled phasegate it is obtained 
from the non-local phase x, 



C = 



■ X 
sin — 
2 



(14) 



In the presentation of our numerical results below we will 
use the motional purity Fbo and the non-local phase x m 
addition to the fidelity F to analyze the performance of 
the optimal pulses. 



III. RESULTS 
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In order to obtain a clear physical picture of the limit- 
ing factors that influence the speed of the two-qubit gate 
operation, we start by exploring a regime in which the 
atom-atom interaction would be so strong to yield a time 
scale shorter than any other in the problem. This regime 
is experimentally unfeasible with optical potential both 
in terms of length scales and of confinement strengths, 
and in this sense the calculation represents little more 
than a toy model. Nevertheless, we simulate the dynam- 
ics taking into account in detail the physical features of 
a real atomic species as if the geometry considered were 
realizable in the laboratory via some trapping force. This 
allows us to gain a thorough understanding of the rele- 
vant energy and time scales, which we subsequently ap- 
ply to a more realistic case of Rydberg-excited atoms 
interacting at longer distances, compatible with realistic 
optical potentials. 



Optimization for Two Calcium Atoms at Short 
Distance 



We consider two calcium atoms in an optical lattice at 
a distance of d = 5 nm that will be excited into a low- 
lying excited state. While such a distance is not feasible 
with the trapping techniques currently available in ex- 
periments, larger distances do not provide a sufficient in- 
teraction strength in the excited state to reach any signif- 
icant fidelity in a reasonable amount of time. Neverthe- 
less, this unrealistic assumption allows us to determine 
the physical mechanisms that limit the gate operation 
time. 

At close distance, the ground state wave functions of 
the two atoms in a harmonic trap can have a significant 
overlap, 



(^ + |*_) « e ^^fr . 

In order to be able to treat qubits carried by the two 
atoms as independent, we compensate for the small value 
of d with an artificially large trap frequency. A choice 
of oj = 400 MHz ensures that the overlap of the wave 
functions is smaller than 10~ 4 . The grid parameters 



FIG. 4: (color online) Fidelity F, non-local phase \ ( m units 
of 7r) and vibrational fidelity, i.e., projection onto the vibra- 
tional target state, Foo, for different gate times T. The inset 
shows the infidelity 1 — F and the respective quantities 1 — x 
and 1 — .Foo The interatomic distance is d = 5nm. 



need to be chosen such that a reasonable number of 
trap eigenstates (about 50 in our case) is correctly repre- 
sented. This is accomplished by taking i? m in to be 5.0 a , 
Rmax = 300.0 ao, the number of grid points Nr = 512, 
with mapping parameters (3 — 0.5, E max = 1 x 10 -8 , cf. 
Ref. dlj. 

The minimum gate operation time to achieve a high- 
fidelity implementation can be due to (i) the strength of 
the exchange interaction in the excited state, or (ii) the 
vibrational motion in the trap. We investigate both hy- 
potheses. The time scale associated with the interaction 
strength is estimated from the maximum phase that can 
accumulate in the interacting state during time T . For a 
non-local phase of one radian, we find 



TZ d (d) = 



1 



Eoa-Voa(d) 



(15) 



where Eo a denotes the energy of two infinitely sepa- 
rated, i.e. non-interacting atoms in the |0a) or |a0) state, 
and Vo a (d) the interaction potential at distance d. For 
d = 5nm, this yields T[£f w 1.23 ps for a non-local phase 
of one radian, and Tf nt « 4.4 ps for a non-local phase of 
7T. The time scale associated with the vibrational mo- 
tion in the trap is estimated by considering the mean 
energy difference of the trap ground state energy to its 
neighboring levels, i.e. the last bound state and the first 
excited trap state. For the chosen trap frequency, we 
obtain T v w 800 ps. 

The optimization results for gate operation times var- 
ied between the two limits T^f and T v are shown in Fig- 
ure |U We compare fidelity F, Eq. ([7]), non-local phase 
X, Eq. flT3J), and vibrational fidelity F 00 , Eq. (JTTJ). The 
optimizations are converged to within AF < 1 x 10~ 4 
except for T = 30 ps and 50 ps which are converged to 
within AF < 2 x 10~ 4 . For durations below 150 ps, with 
errors remaining larger than 10~ 2 no satisfactory fidelity 
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is obtained. As the gate operation time approaches T v , 
optimization is successful in the sense that fidelities arbi- 
trarily close to one can be reached. The results shown in 
Fig-Hlcan be understood as follows: The two-qubit phase 
X increases with the pulse duration T, and at T = 5 ps, 
the time that was roughly estimated to reach a non-local 
phase of it, about half that phase is actually obtained. 
This is not surprising since the wavepacket is not in the 
excited state for the complete gate duration T due to 
the switch-on and switch-off phases of the pulse and its 
general shape. The non-local phase reaches the desired 
value of 7r at about 50 ps. We thus find that a prolonged 
action of the exchange interaction leads indeed to a non- 
local gate. However, for short gate durations, no control 
over the motional degree of freedom can be exerted. For 
T = 5ps, the vibrational fidelity drops below 50%, and 
it increases rather slowly for larger T . This is due to the 
wave packet spending enough time in the excited state 
to be accelerated by the 1 / B? potential. When the laser 
pulse returns the wavepacket to the electronic ground 
state, it has acquired significant vibrational energy. Since 
the pulse is too short to resolve the vibrational motion 
in the trap, optimization cannot identify the desired trap 
state and thus it cannot counteract the excitation. Popu- 
lation of excited trap states after the gate can be avoided 
only once the pulse is long enough to resolve different 
trap states. As the gate duration becomes comparable 
to T v , fidelities close to one are obtained. 

An analysis of the dynamics induced by the optimized 
pulses is instructive for short gate durations despite the 
low fidelity. Figures [5] and [5] display the optimized pulse, 
its spectrum, and the dynamics for T = 5ps. The 
guess pulse that is used to start the iterative optimiza- 
tion has a Gaussian envelope with a peak intensity of 
about 4.9 x 10 7 V/m. The intensity was chosen to drive 
one complete Rabi cycle for a single qubit in the |0) state 
(2-7r-pulse). The pulse fluence is increased by a factor of 
about 7 during the course of iterations. Optimization re- 
sults in a pulse shape that clearly shows more features 
than a Gaussian, cf. top panel of Fig. [SJ The first peak 
of the pulse, centered around w 100 fs, drives significant 
population transfer to the auxiliary state. The last peak 
of the pulse, centered around w 1 ps, restores all pop- 
ulation to the electronic ground state. The dynamics 
yielding the non-local phase occur at intermediate times. 
Since there is no interaction between the two atoms in 
the |01) and \al) states, the population dynamics of the 
|01) state (red dot-dashed line in Fig. [5]) are equivalent 
to those of a single atom. Comparison of the popula- 
tion dynamics of the |00) state and the |01) state (black 
solid and red dot-dashed lines in the top panel of Fig. [5]) 
therefore yields the difference between single-qubit and 
two-qubit dynamics. For the short gate operation time 
shown in Fig. [5j the two curves are fairly similar. This 
is in agreement with the non-local phase of only 0.63 it 
that is achieved. The spectrum of the optimized pulse, 
cf. bottom panel of Fig.[5J is tightly centered around the 
|0) — > | a) transition frequency of 23652.30 cm" 1 . It is 
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FIG. 5: (color online) Optimized pulse (gray) and population 
dynamics ( 1 00) state: solid black line, |01) state: dot-dashed 
red line, top) and pulse spectrum (bottom) for T = 5 ps (F = 
0.805). 



sufficiently narrow to guarantee that no undesired transi- 
tions, for example into the 1 H g (S + D) state, are induced. 

The dynamics of the two-qubit system and of a sin- 
gle qubit can be analyzed by projecting the time-evolved 
two-qubit basis state onto the initial two-qubit and 
single-qubit states, 

Tij (t) = (y(Je)|0(t,0;e**)|y(iJ)) , (16) 



U(t,0;e opt ) 



(17) 



The phase dynamics Tjj(i) and Tj(t) obtained with the 
optimized pulse of Fig.[5]are shown in Fig.[6]with the top 
(bottom) corresponding to the two-qubit (single-qubit) 
dynamics. The phase of the initial state is indicated by a 
filled black circle, the phase of the final state by a black 
square, and the phase of the optimization target by a 
blank circle. The phase dynamics of the j 01} state (up- 
per right panel of Fig. [SJ no connecting lines are shown) 
consist of a mixture of the natural time evolution and the 
dynamics induced by the pulse. The |1) state (lower right 
panel of Fig. ^ displays only the natural time evolution. 
Since it is not coupled to any other state by the pulse, 
the population remains one, i.e., on the unit circle, at all 
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FIG. 6: Phase dynamics induced by the optimized pulse 
(T = 5 ps) in the complex plane for the two-qubit states (top) 
and the single-qubit states (bottom). The initial state is rep- 
resented by a black circle, the final state by a black square, 
and the target state by a blank circle. 



times. The fact that the 1 00) state does not return to the 
unit circle at the final time T indicates leakage out of the 
quantum register. The overall fidelity amounts to only 
0.805, cf. Figure El with the target phase for both the 
|00) state and |01) state missed by almost equal amounts, 
cf. black squares and empty circles in the upper panel 
of Fig. [5] This reflects that the optimization is balanced 
with respect to all targets, i.e. the terms in the sum 
of the target functional, Eqs. ©-(O, all enter with the 
same weight. A comparison of the |00) and |0) phase dy- 
namics (upper and lower left panels of Fig. [5]) illustrates 
how a true non-local phase is achieved, even though the 
optimization is only partially successful: Without inter- 
action the phase on the |00) state would evolve according 
to 0oo — 20o- The extent to which this is not the case 
demonstrates how the interaction leads to the non-local 
phase. 

The optimized pulses, their spectra, and the corre- 
sponding population dynamics for intermediate and long 
gate durations are shown in Figs. [7] (T = 50 ps) and [5] 
(T = 800 ps). The guess pulses were again chosen to 
be Gaussian 27r-pulses. During the course of the itera- 
tions, the pulse fluence was increased by a factor of 28 
for T = 50 ps and by a factor of 44 for T = 800 ps. The 
overall structure of the optimized pulse for T = 50 ps is 
similar to that obtained for T = 5ps: Two peaks at the 
beginning and the end induce population transfer to and 
from the auxiliary state while the intermediate part of 
the pulse drives Rabi oscillations in the course of which 
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FIG. 7: (color online) Pulse dynamics (top) and spectrum 
(bottom) for the optimized pulse with T = 50 ps after 255 
iterations(F = 0.988), analogously to Fig. [5] 



the non-local phase is achieved. The first peak of the 
pulse triggering population transfer to the auxiliary state 
remains clearly visible as the gate operation time T is fur- 
ther increased, cf. Fig. [8] Overall, however, the optimal 
pulse shows less discernible features for T = 800 ps than 
for the shorter gate durations where a sequence of sub- 
pulses was found. This is reflected in the population dy- 
namics: While for T — 50 ps, each subpulse drives a par- 
tial transfer, resulting in step-wise population dynamics, 
an almost adiabatic behavior is observed for T = 800 ps. 
Comparing the population dynamics for a single qubit 
and the two-qubit system for T = 50 ps (black solid and 
red dot-dashed lines in the upper panel of Fig. [7]), more 
differences are obtained than for T = 5ps, cf. Fig.[5j but 
overall the single qubit and two-qubit dynamics are still 
fairly similar. This changes dramatically for T = 800 ps 
(black solid and red dot-dashed lines in the upper panel 
of Fig. [5$, where the populations dynamics for 1 00) and 
j 01) are clearly distinct, reflecting that the desired non- 
local phase is fully achieved (x = 0.9987T for T = 800 ps 
as compared to x = 0.9757T for T = 50 ps). The spectrum 
of the optimal pulse for T = 50 ps is fairly similar to that 
obtained for T = 5ps, cf. the lower panels of Figs. [S] and 
[7] It basically consists of a single narrow peak centered 
around the |0) — > \a) transition frequency. The spectrum 
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FIG. 8: (color online) Pulse dynamics (top) and spectrum 
(bottom) for the optimized pulse with T = 800 ps after 104 
iterations (F = 0.999), analogously to Fig. [5] 



for T = 50 ps shows somewhat more features within the 
peak which is attributed to the better spectral resolution 
for larger T. For T — 800 ps, the spectrum of the op- 
timal pulse consists of a narrow peak at the |0) — > \a) 
transition frequency and sidebands. These sidebands re- 
main sufficiently close to the |0) —> \a) transition, that 
resonant excitation into other electronic states can be ex- 
cluded, cf. Fig. [3] The phase dynamics induced by the 
optimized pulse of Fig. [5] are shown in Fig. |H1 Since the 
target phases include the natural time evolution, their lo- 
cations in Fig. differ from those in Fig. [5] The overlap 
of the final states (black square) and the target states 
(open circle) confirms success of the optimization. All 
phases end up on the unit circle demonstrating that no 
leakage from the quantum register occurs at the end of 
the gate for T = 800 ps. 

We also carried out optimizations for non-local target 
phases that are a fraction of it such as | or |. If high- 
fidelity implementations of such fractional phasegates are 
found, several of these gates can be combined sequentially 
to yield a total non-local phase of n. However, for short 
gate operation times, optimization for non-local target 
phases smaller than tt did not prove any more success- 
ful than optimization for w. In particular, population of 
excited trap states at the end of the gate could not be 
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FIG. 9: Phase dynamics induced by the optimized pulse (T = 
800 ps) in the complex plane for the two-qubit single-qubit 
states, analogously to Fig. [6] 



avoided for fractional phasegates either. Moreover, we 
investigated whether pulses driving multi-photon transi- 
tions, for example pulses with their central frequency a 
third of |0) — > |1) transition frequency, yield better fi- 
delities for short gate operation times. However, we did 
not observe any substantial difference in the results com- 
pared to the pulses reported in Figs. [5l [7j and [8] These 
additional investigations confirm that for our example of 
two ultracold calcium atoms in an optical lattice, the lim- 
its on the gate operation time is set by the requirement 
to restore the ground vibrational state of the trap. 



B. Optimization for Two Atoms at Long Distance 
under Strong Dipole-Dipole Interaction 

To determine whether it is really the ground state mo- 
tion in the trap and not the non-local interaction in the 
excited state that sets the speed limit for two atoms res- 
onantly excited to an interacting state, we vary the in- 
teraction strength C3 of the dipole-dipole interaction po- 
tential, 



V(R) 0a = V(R) a0 = -^, 



(18) 



keeping the trap frequency constant. We consider the 
atoms to be separated by d = 200 nm which corresponds 
to a realistic optical lattice in the UV regime. In order 
to keep the overlap of the ground state wave functions 
smaller than 10~ 4 at a distance of 200 nm the trap fre- 
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converged to AF < 1 x 10 -4 The interatomic distance is 
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quency has to be set to at least 250 kHz. This corre- 
sponds to T v rj 2 ns. 

For the interaction potential of two calcium atoms in 
the B 1 ^^ state used in Sec. lIII Al the C3 coefficient takes 
a value of 16.04 a.u.= 0.5217 x 10 3 nm 3 cm- 1 0, [H, [H| . 
This results in an interaction energy of about 4 cm" 1 at 
d = 5 nm. Based on the results of Sec. MI Al we know 
that such an interaction energy is sufficient to yield a non- 
local phase in a few tens of picoseconds. For d = 200 nm, 
the same interaction energy is obtained by choosing C 3 
to be roughly 1 x 10 6 a.u. Just for comparison, the C3 
coefficient for highly excited Rydberg states is about 3 x 
10 6 a.u., resulting in an interaction energy of about 1.3 x 
10 _3 cm _1 at a typical distance of 4/xm for two atoms 
trapped in optical tweezers [32j. 

We vary the C3 coefficient from 1 x 10 6 a.u. to 
1 x 10 9 a.u. If the gate duration is solely determined 
by the requirement of a sufficiently strong interaction 
to realize the non-local phase, we expect to find high- 
fidelity implementations with optimal control by increas- 
ing the C3 coefficient. In particular, we pose the question 
whether picosecond and sub-picosecond gate durations 
can be achieved given that the interaction is sufficiently 
strong, i.e. given that the C3 coefficient is sufficiently 
large. Based on the results of Sec. IIII Al where a non- 
local phase of n was achieved within 50 ps, we estimate 
that C3 needs to be increased from 1 x 10 6 by a factor of 
50 (100) to obtain a high-fidelity gate for a duration of 
lps (0.5 ps). 

Figure [TUl presents optimization results for a controlled 
phasegate with gate operation times of T = 0.5 ps and 
T = 1 ps. The central frequency of the guess pulse 
was adjusted in each case to compensate for the in- 
creased interaction energy and ensure resonant excita- 
tion. The grid parameters were chosen to be i? m i n = 5 ao, 
i? max = 13000 a w 688 nm, and N R = 2048. This choice 



of i? m ax guarantees that at least fifty eigenstates of the 
trap are accurately represented. We verified that the grid 
is sufficiently large, i.e., the wave packet does not reach 
the boundaries of the grid during propagation. Moreover 
we checked that doubling the number of grid points did 
not yield substantially different results. Figure [TUl clearly 
shows that increasing the interaction strength leads to 
larger non-local phases. A non-local phase of ir is reached 
for C3 = 4 x 10 8 a.u. However, increasing the interaction 
strength also results in a complete loss of control over 
the motional degree of freedom, with the vibrational fi- 
delity F00 reduced to 75% and below. A stronger inter- 
action in the excited state accelerates the wave packet 
more, increasing its vibrational excitation. This results 
in a larger spread over many trap states upon the wave 
packet's return to the ground state. Since 0.5 ps or lps 
are much, much shorter than the time scale of the mo- 
tion in the trap, T v , the optimization algorithm cannot 
resolve the eigenstates of the trap. Thus it cannot iden- 
tify the target Fbo which is consequently missed com- 
pletely. Note that trap frequencies up to a few MHz are 
possible and imply shorter T v . Nevertheless, any realis- 
tic trap frequency results in vibrational time scales much 
larger than a few picoseconds. While excited state po- 
tentials providing for a strong interaction between two 
neutral atoms exist, resonant excitation into such an ex- 
cited state will not yield an ultrafast non-local gate due 
to coupling with the motional degree of freedom, unless 
we consider a gate scheme that is completely insensitive 
to vibrational excitations 



IV. SUMMARY AND CONCLUSIONS 

We have studied high-fidelity implementations of a 
controlled phasegate for two trapped ultracold atoms via 
resonant optical transitions to an electronically excited 
state with long-range diatomic interaction. To the best 
of our knowledge, we have for the first time explicitly ac- 
counted for the detailed independence of the interaction 
and thus for the coupling between electronic and nuclear 
dynamics that may cause leakage out of the quantum 
register. We have employed optimal control theory to 
calculate laser pulses that carry out the gate. This has 
allowed us to determine gate implementations of basically 
arbitrarily high fidelity provided the gate operation time 
is sufficiently long (and at the same time short enough 
to neglect dissipation). Our main goal was to achieve 
the fastest possible gate implementation and to identify 
what limits the gate operation time, i.e., to determine 
the quantum speed limit for a controlled phasegate for 
two neutral trapped atoms. 

The standard reasoning considers the interaction 
strength to be the limiting factor, i.e., the gate opera- 
tion time is estimated by the inverse of the two-qubit 
interaction. Our calculations show that a second time 
scale might come into play: For resonant excitation, the 
interaction between the two atoms causes a coupling be- 
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tween electronic and nuclear dynamics. This induces vi- 
brational excitation which can be carried away by the 
laser pulse only if the target state is fully resolved during 
the optimization. The gate operation time is thus lim- 
ited cither by the two-qubit interaction strength or by 
the vibrational motion in the trap, which ever one of the 
two yields the larger time. 

This finding has important implications for the design 
of two-qubit gates where the qubits are carried by neu- 
tral atoms. For example, excitation of atoms into Ryd- 
berg states yields an interaction that one might expect to 
allow for nanosecond to sub-nanosecond gate operation 
times. However, the motional state of the atoms needs 
to be restored at the end of the gate, and traps with sub- 
nanosecond vibrational motion seem difficult to realize. 
Thus the question of how an ultrafast two-qubit gate can 
be realized in a scalable setup still remains open. 
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Appendix A: Reduced Optimization Scheme 

Since the dynamics relevant for obtaining a non-local 
phase involves only the |00) state out of the four two- 
qubit states, cf. Fig. [2] we can reduce our full model, 
Eq. to one describing only the left-most column of 
Fig. [2j This is a direct consequence of optimizing for 
a diagonal two-qubit gate. We are then operating in a 
4iV#-dimensional Hilbert space instead of a 3 x 3 x Nr- 
dimensional Hilbert space. However, care must be taken 
to extract the correct non-local phase x- The phase 0oo 
describing the time evolution of the two-qubit |00) state 
alone is not sufficient to obtain the non-local phase X 



since 0oo contains contributions from both the single- 
qubit and two-qubit dynamics. We therefore need to aug- 
ment our reduced model for the dynamics starting from 
the |00(iZ)) by a two-level system (|0), \a)). This captures 
the purely single-qubit dynamics in the phase 0o- The 
non-local phase is then obtained as the difference between 
0oo and (twice) 4>o, see Table HI The optimization targets 
for the full model and the reduced model are also listed in 
Table HI They correspond to the unitary transformation 
for the full model and to two state-to-state transitions 
for the reduced model. Note that this type of state-to- 
statc transition requires a phase-sensitive functional such 
as the one in Eq. (JT]). We have checked numerically that 
the full and reduced model are indeed equivalent: Prop- 
agating the Schrddinger equation with an optimal pulse 
obtained for the reduced model but employing the full 





full reduced 


target 


|00) ->• e'^+M | 00 ) 

|01) -> jOl) |00} -> e l W+<« |00> 
|10) ->e 10T |10) |0) -> e 4 " 1 -' 2 |0) 
|11) -> e l4>T |11) 


gate phases 


000 = 000 

4>w = 4>oi = <f>o + <t>i 
0n = 2<t>i 


non-local 
phase 


X — 0oo — 0oi — 0io + 0ii X = 0oo — 20o 



TABLE I: Comparison between the full model and the re- 
duced model for optimization of a controlled phasegate. 



Hamiltonian, Eq. ©, we obtained the same fidelity as 
for the reduced model. 
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